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Abstract 

The present work revisits the problem of modeling diffusion above a stagnant liquid interface in a 
capillary tube geometry. In this revisitation we elucidate a misconception found in the classical model 
proposed by Bird et. al. Furthermore, we propose alternative explanations for thermally forced diffusion 
and provide a description of natural convection in the absence of forcing terms. 

1 Introduction 

The Stefan diffusion tube problem is a well studied phenomenon that models the diffusion of a species 
through an ideal binary gas mixture above a hquid-gas interface [2j|4j[6|[l0j[T3j[28|[29] . Researchers typically 
couple Fick's first l&w with a mass balance equation to arrive at a transport equation for the evaporating 
(or condensing) species (eg water vapor). A classical model was presented by Bird, Lightfoot, and Stewart 
in 1960 (second edition published in 2007) [4] and involves several assumptions regarding the density of the 
bulk gas and the relative fluxes of the species. In the present work we elucidate the assumptions made and 
discuss a possible misconception with this classical model. 

The Stefan diffusion tube problem can be extended to describe diffusion through a porous medium. It is 
well documented that diffusion through porous media can act faster than that predicted strictly by Fick's first 
law (8)[9)[18)[22}|27|. In 1957, Philip and deVries (8)[9)[l8] postulated a description of enhanced vapor diffusion 
through porous media via gradients in temperature and saturation. These gradients do not naturally occur 
with the classical combination of Fick's first law and mass balance, and this model has come under recent 
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scrutiny 21 22 24 25 27 due, in part, to the difficulty in measuring pore-scale effects. In the present work 
we suggest an alternative explanation to the effects of thermal gradients on vapor diffusion. 

In Section [2] we give a thorough derivation of the mass transport equation and propose several alternative 
dependent variables (mass concentration, relative humidity, and chemical potential) for re-expressing Pick's 
law. The resulting models take the form of advection-diffusion equations, so the Peclet number naturally 
arises as a dimensionless coefficient. In Section [s] we focus of diffusion dominated systems (Peclet number 
much less than 1). 

As part of the discussion of diffusion dominated systems in Section [3] we discuss the model proposed by 
Bird, Lightfoot, and Stewart [4j. They formulated a nonlinear equation for steady state diffusion in a binary 
gas mixture. In their derivation it was noted that the velocity (and hence the flux) of the inert species 
in the mixture was negligible as compared to that of the diffusing species. In section [3. 1| we will re-derive 
the Bird model in terms of a transient diffusion equation, but here we will not discard the velocity of the 
inert species. It will be shown that a simple momentum balance can be used to revert the classical model 
back to a linear diffusion model. Section |3.2| discusses the necesesity of using a time dependent equation for 
certain length scales versus using a steady-state (or quasi steady-state) model as proposed in j4| . In Section 



3.3 we propose another model for thermally enhanced vapor diffusion as compared to [8j[9j[18]. This model 
verifies that pore-scale gradients in temperature can cause thermally enhanced diffusion (as predicted), but 
the proposed mechanism is different than that of the classical models. 

In Section|4]we seek to couple the mass and momentum equations to form a system of governing equations 
where the Peclet number is non- negligible. Whitaker 29 1 proposed using momentum balance equations to 
describe the flux of the species in a Stefan diffusion tube, but in the present work we pair the mass and 
momentum equations to form a coupled system of governing equations. 

2 Derivation of Advection Diffusion Models 

The standard advection diffusion equation is a combination of Pick's Law for diffusion and the mass balance 
equation. Por a system consisting of an ideal mixture of water vapor, g^^ and inert air, ga, Pick's Law can 
be written as 

JS" = -DVpf", (2.1) 

where J^" is the mass flux of water vapor [kg(gt,)/(m^- s)], D is the diffusivitiy [m^/s], and p^" is the mass 
density (concentration) of water vapor \kg{gy) /w?]. For a complete notation reference see Appendix [Aj The 
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mass balance equation for species j in the gas phase can be written as 

+ V • (p«vs^ ) = fs^- (2.2) 

dt 

where v^j is the velocity of species j within the gas mixture relative to a fixed frame of reference, and f^^ is 
a mass exchange term accounting for chemical reactions between species [29] . In the present work we assume 
that — 0. The combination of these two equations for water vapor (j = v) gives a transport equation 
for the mass of water vapor via advection and diffusion. The remainder of this section gives details of the 
derivation, a discussion of the choices of dependent variables, and addresses some possible issues. 

The mass flux, J^" , is a measure of the diffusive flux of the water vapor through the mixture. For that 
reason, we note that the flux of water vapor can be written as 

Jff" = p9- (v9" - V9) EE p9" vS"'9. (2.3) 

The immediate concern when coupling (??) with (??) is that the mass balance equation is written in terms 
of the absolute velocity of the water vapor, v^" , and the flux is written in terms of the diffusive velocity, 
v^"'^. Typically [3||4]|29], the bulk velocity of a mixture is defined as a mass (or mole) weighted. Here we 
use a mass weighted velocity so that the flux of the bulk gas is expressed as the weighted sum of the flux of 
each species, p^v^ = SjLi P^^v^^ . In a binary mixture this takes the form 

p9v9 =p9"v9" +pf"vS». (2.4) 

The mass flux of the water vapor can be written in terms of the diffusive flux by adding and subtracting 
the bulk flux to arrive at 

p9"vf" = pS-vf-f + pf"vf = Jf" + p9"v9. (2.5) 

Substituting this into the mass balance equation, (??), writing the flux in terms of Pick's Law, (??), and 
rearranging gives an advection-diffusion equation with non-constant velocity and diffusivity. 

^ + V-(pS^'v9) = V-(DVp»"). (2.6) 

We have not made the assumption that D is constant in space. While it may be true that D is constant in 
space in many application it is worth noting that Z? is a function of temperature. Therefore, when coupling 
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with thermal effects, the diffusion coefficient must be left inside the divergence operator. See Section 3.3 for 
a discussion on thermal effects and non-constant diffusivity. 



2.1 Relative Humidity Model 

To arrive at a model in terms of relative humidity we begin by nondimensionalizing equation (??). In this 
case we assume that D is constant in space so it can be factored from the right-hand side (we will return to the 



temperature dependence in Section 3.3). To nondimcnsionalize we let L be the length of the capillary tube 
(the diffusive length of interest), and scale the spatial and temporal variables as x — L^^ and t — [LF' /D)t 
respectively. To give a dimensionless form of (??) we need to scale the density, p^" , and the bulk velocity, v^, 
via some characteristic values. Presuming that the temperature is constant throughout the diffusive domain, 
we let p^" — psat{T)ip, where Psat{T) is the saturated vapor density at the given temperature. Observe 
that this definition (along with the ideal gas law) defines ip as the relative humidity (divided by 100). For 
the bulk velocity we scale by a problem dependent characteristic velocity so that = WcV^. The value 
of Vc depends on the experiment of interest. Under the assumption that the temperature is constant, the 
advection diffusion equation is written in terms of relative humidity as 

^ + (Pe)Ve • (^v^) = • Ve^. (2.7) 

This is only one choice of dimensionless parameters, but the relative humidity is a natural choice as it is 
easy to measure experimentally. 

The coefficient of the advection term is the Peclet number and is formally the ratio of the advective 
velocity to the diffusive velocity: 

Pe= — - = , , (2.8) 
D (D/L) ^ ' 

If a system is diffusion dominated then Pe <^ 1 and this term can possibly be neglected. On the other 
hand, if there is some driving force acting on the bulk velocity then Pe > 1 and the advective term is 
non-negligible. A similar model in terms of relative humidity (without the advective term) has been used by 



D. Or, E. Shahraeeni, and colleagues for pore scale modeling of evaporation and condensation 20 -22 . 



2.2 Chemical Potential Model 

In thermodynamics and physical chemistry [5| |15| the chemical potential is the natural descriptor of mass 
diffusion. Matter is known to move from areas of high chemical potential to areas of low chemical potential, 
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and at equilibrium the chemical potential of a species in the gas phase is equal to the chemical potential of 
that species in the liquid phase [l][5{ [l2p^ . These descriptors of the chemical potential make it an appealing 
dependent variable for modeling diffusion dominated systems. 

Using the ideal gas law [5 15 , the chemical potential is defined in terms of the partial pressure of the 



water vapor in the gas phase and a reference chemical potential; 

=f,3^ +R3^T\n(^^^ (2.9) 

(a subscript * represent reference quantities, see Appendix[A|). Solving for the ratio of pressures and recalling 
that the relative humidity is defined a,s ip ^ P^"' /Psat we arrive at a conversion between relative humidity 
and chemical potential in the gas phase: 

e" = Xifi, (2.10) 

where u = (fi^^ — /if^')/(i?^''T) and A = Psat/p*- From equation (??), Pick's law can now be rewritten as 

39^ = -DVpS^' = -DV iPsat{T)ip) = - (^^^Pf^^ Ve". (2.11) 

Using the definition of mass weighted velocity (??), the advection-diffusion equation for relative humidity, 
(??), and cancelling like terms gives an advection-diffusion type equation in terms of the dimensionless 
chemical potential 

de" 

— + (Pe)V5 • (e"v9) = • V^e". (2.12) 

OT 

Setting the Peclet number to zero (neglecting advection) gives a diffusion equation in terms of the dimen- 
sionless chemical potential. Curiously, though, this equation does not directly solve for chemical potential. 
Of course, from a practical standpoint there is no problem solving u. The point, though, is that the chem- 
ical potential is no better mathematically suited to describe diffusion in a gas as relative humidity or mass 
density. 

Since each of the three proposed models for diffusion of a gas ((??), (??), and (??)) are roughly the same 
mathematically, it remains to determine (a) how these models compare with classical diffusion models such 
as that of Bird, Lightfoot, and Stewart, (b) whether the transient nature of the equations are relevent for 
certain length scales, (c) the sensitivity of the diffusion to thermal gradients, and (d) the role of the advective 
term. The next sections address these issues. For simplicity we will mainly focus on the relative humidity 
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form of the advcction diffusion equation, (??). 



3 Diffusion Dominated Models 

In the Stefan diffusion tube problem the physical action is assumed to be dijfusion dominated. That is to 
say that the diffusion is the dominant physics and advcction is either absent or negligible. A classical model 
for this problem was presented by Bird, Lightfoot, and Stewart in 4| and has subsequently been used by 
several other authors. This section discusses the Bird et. al. model along with relevant temporal scales and 
thermal effects for diffusion dominated systems. 

In a diffusion dominated system the Peclet number is small, and hence the advective term, (Pe) • {(fv^), 
is negligibly small as compared to the time derivative and the diffusion term. Under this assumption, (??) 
becomes the standard diffusion (or heat) equation. 



d(p 



(3.1) 



Consider now a capillary tube geometry with a stationary liquid-gas interface and experimentally controlled 
ambient conditions (see Figure [TJ. Assume that the interaction between the gas phase and the solid walls 



9v + 9a 



Figure 1: Cartoon geometry of a typical capillary tube. At ^ = the Kelvin Equation (??) gives the 
boundary condition. At ^ = 1 the boundary condition is assumed to be under experimental control. 



of the tube is negligible so that the diffusion is predominantly one dimensional up the vertical axis of the 
tube. This assumption is confirmed to be approximately valid for water vapor diffusing through a tube 



by Salcedo-Daz, Ruiz-Femenia, Kerkhof, and Peters 19 . Take L to be the length of the tube and define 



^ = to be the liquid-gas interface (hence defining the geometry as ^ € [0, 1]). Given a tube diameter, the 



liquid-gas interface boundary condition is given by the Kelvin equation 16 , 17 20 



<Plc=o = exp 



— 2(TK 



(3.2) 



Here, a is the surface tension [N/m], k is the mean curvature of the interface [l/m], f?^" is the specific gas 
constant for water vapor [J/(kg • K)], and p' is the hquid density [kg/m'^]. This boundary condition expresses 
the balance of forces at the hquid-gas interface. Given a one dimensional mathematical geometry we can 
see the capillary tube as having a flat interface with the boundary condition governed by curvature and 
surface tension parameters depending on the actual three dimensional experimental setup. The boundary 
condition on the open end of the capillary tube could be either Dirichlet-type (fixed relative humidity) or 
Neumann-type (fixed fiux). For either set of boundary conditions (??) has a standard Fourier series solution. 

In the following subsections we comare with the Bird, Lightfoot, and Stewart model, discuss relavent 
time scales for the temporal solution, and discuss the sensitivity to thermal gradients. Throughout these 
subsections we compare with solutions to the linear diffusion equation (??). 

3.1 A Comparison with Bird, Lightfoot, and Stewart 

In [4], Bird, Lightfoot, and Stewart examined the same problem of the Stefan diffusion tube. Their approach 
to modeling this problem (rewritten in the present notation) is to first note that the bulk velocity is mass 
weighted, p^v^ = J2j=v a and to use this to rewrite the mass flux of the water vapor: 



Here we have used the notation C^^ — p^^ / p^ [-] to represent the mass concentration of species j in the gas 
phase. Note here that the mass flux has been written in terms of the diffusive flux and the velocity of the 
dry air relative to a fixed frame of reference. This is in contrast to the derivation of equation (??) where we 
write the mass flux in terms of the diffusive flux and the bulk velocity. 

Using Fick's Law, (??), for p9^\-9^^a^ a^d eliminating pf"v^" in (??) with (??) gives 



pav^av.a pa. {c^.^a. _^ C<^^^a.) . 



Solving for pS^'-va. gives a description for the mass flux of the water vapor: 




(3.3) 




(3.4) 
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Converting the density terms to relative humidity (assuming that temperature is constant) gives 



dip 



+ V • (</9V»<") = V • 



D 



V(/3 



(3.5) 



To simpHfy the fraction 1/(1 - C^-) we first recall that C^" = p^" / and C^- + C^- = 1. Therefore, 



1 pB p3- + p3- 



o9^ 



Psatf 



1 - Cf- C9» p9 



= 1 + — =1 . 

p9a pg. 



At this point we must simpliiy p^<» using the ideal gas law. We make the further assumption that is fixed. 
Dalton's law states that the bulk pressure is the sum of the partial pressures in an ideal gas mixture. If we 
assume that the bulk pressure is fixed then 



(3.6) 



Next use the ideal gas law to rewrite p^^ = R^^ Tp^' ; 



p, = R3-Tp3^ +R3"Tp,atV'- 



(3.7) 



This implies that 



P* R^^Psat 



Ra-T 

\R9aTJ 



-if: 



P* 

f Psat' 

\ p* 



Psat 



Rs^T' 



(3.8) 



Now we see that 1/(1 — C^") can be written as 



1 



1 + 

1 + 
1 + 
1 + 



RP^2^\ f if 

P* J \'^->^<P 
RP^R?^Tp^\ / y 

R^^P* ) \l-Xip 

R^^Psat \ f y 
R3-\\ / if 

R9- J \l- Xtp 



(3.9) 



Rewriting (??) with this simplification and nondimensionalizing finally results in a nonlinear advection 
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diffusion equation (as compared with the linear advection diffusion equation (??)) 

(3.10) 



g + (Pe)V5.(^v«i=V5 



where A —Psat/p* and A ~ (Ai?f°)/i?3" are dimensionless groups. 

It should be noted here that equation (??), with the advective term set to zero, is the transient version 



of Bird's equation |4 . Whitaker 29 suggested that "one can develop convincing arguments in favor of 
..." neglecting the air-species flux term. In the present notation this amounts to arguing that v^" is 
approximately zero (or that Pe <C 1). If we allow then we have fixed a frame of reference such 

that the dry air is viewed as immobile and the water vapor is allowed to diffuse through this stagnant air. 
The trouble here is that this frame of reference is unachievable by an experimentalist. On the other hand, 
a frame of reference where the bulk velocity, v^, is fixed is achievable. 
Given that the weighted sum of the diffusive velocities is zero, 

C-s^v^-s + C9»v»-f = 0, (3.11) 

we know that as the water vapor molecules move upward the air molecules must move downward. This 
indicates that if we were to take the frame of reference where v^" « then the observer would be moving 
downward with the dry air (in an evaporative process) while the water vapor moves upward. 

The Bird, Lightfoot, and Stewart model gives a solution to the diffusion process in this frame of reference, 
but as mentioned prior, the experimentalist cannot achieve this frame of reference. A sample solution to (??) 
(with Pe ~ 0) is shown in Figure [2] (blue) with a comparison to linear diffusion (red); each as viewed from 
the frame of reference where « 0. For this sample solution we considered Dirichlet boundary conditions 
(pi) ^ 1 and iy9i = to best show the perceived enhancement. A jump-type initial condition was used with 
(^(^ = 0,r = 0) = (^0 and </j(^ 7^ 0,r = 0) = 0. MATLAB's pdepe routine was used to solve the nonlinear 
model. In the enlarged portion of Figure [2] it is clear that the nonlinear diffusion model is advancing faster 
in time than that predicted by the linear model. In Figure [3] we see that the steady state solutions for the 
two equations differ by approximately 10~^. 

Our conclusion for this subsection is that it is possible to mis-interpret the results of the Bird et. al. 
model. In order to properly use this model while viewing an experimental apparatus one must not neglect 
the advective term. This introduces a closure issue to the model as v^" is now an unknown. In contrast 
with Whitaker [29[ , we choose to leave the advective term in the equation and rewrite using equation (??) 
as a simple expression of momentum balance. With (??) and Fick's Law we can rewrite the advective term 
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Time Evolutions of Linear & Nonlinear Diffusion IVIodeis 




Figure 2: Perceived enhancement in diffusion with the model of Bird et. al. 4J. Observe that the nonlinear 
model (blue) is evolving faster in time than the linear diffusion model (red). Simulation: Tmax — 0.25 with 
tolerance 10"® for MATLAB pdepe. 

in terms of and Vp^" as follows: 



D 



1 \ / D 



ps J V 1 - Cf" 



(3.12) 



Replacing the v^" term in (??) yields 



Qp9. ( (79^. 



dt 



1 - 



1 



1 - Cf" 



After simplifying we return to the linear advection diffusion equation (??) with D held constant in space 



9pf 



dt 



( 1 - Cf- 
V • (p^"v^) = W • (^^-^Vp^^ 

= L'V • Vp9" . 



(3.13) 



No additional assumptions went into this derivation, and therefore, this provides another (admittedly con- 
voluted) derivation for the linear advection diffusion equation. It also elucidates the fact that the nonlinear 
diffusion equation (??) falsely reports enhanced diffusion as compared to the linear diffusion model. This 
point was not made in [4], and the advective term was neglected leading to the nonlinear diffusion model 
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Time Evolutions of Linear & Nonlinear Diffusion Models 



Comparison on Linear and Nonlinear Diffusion 




Nonlinear vs. Linear 
Linear vs. Linear Steady State 
Nonlinear vs. Nonlinear Steady State 
Nonlinear vs. Linear Steady State 



0.4 0.6 




Figure 3: Perceived enhancement in diffusion with the model of Bird et. ah [4]. Left: transient solution 
of both diffusion models. Right: relative errors. Observe that the relative error between the nonlinear and 
linear models is C(10"^). Simulation: T^^ax — 1 with tolerance 10~^ for MATLAB pdepe. 



(nondimensionalized and put into the present notation) 



(3.14) 



Bird's model has been used in several works without regard to the momentum balance proposed here. 
Examples of such works include (6)[TlJ[l3)[l4l[l9)|28}|3l] . 



3.2 Relevant Temporal Scales 

We return now to the linear diffusion dominated system modeled in terms of relative humidity; equation 
(??) with Pe <C 1. Our focus in this subsection is to determine when a steady state vs. a transient solution 
is necessary to model diffusion. The fact that (??) is in dimensionless form hides the effects of the scale of 
the problem (L) on the transient nature of the solution. Table [I] shows the time needed to reach within 1% 
L2 relative error of steady state given various lengths of capillary tubes: 



E(r) 



\\y^ii,T) - (p{^,T = +00) 

\\ipi^,T = +00)||2 



(3.15) 



Notice that for longer tubes the transient nature of the solution plays a role for longer time and therefore 
cannot be neglected. For shorter tubes, on the other hand, the transient solution approaches the steady 
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(sec) 


(sec) 


(sec) 






L — Im 


L = 1cm 


L ~ Imm 




0.44 


1.71 X 10"* 


1.71 X 10" 


1.71 X 10-2 


ipi = 0.2 


0.41 


1.60 X 10^ 


1.60 X 10" 


1.60 X 10-2 


ipi = 0.4 


0.37 


1.44 X 10-* 


1.44 X 10" 


1.44 X 10-2 


ipi = 0.6 


0.32 


1.25 X 10^ 


1.25 X 10" 


1.25 X 10-2 


ipi = 0.8 


0.24 


9.34 X 10^ 


9.34 X 10-1 


9.34 X 10-3 


<pi = 0.95 


0.09 


3.50 X 10^ 


3.50 X 10-1 


3.50 X 10-3 



Table 1: Time to reach within 1% L2 relative error of steady state for a diffusion dominated system with 
capillary tube length L. Values are found by setting equation (??) equal to 0.01 and numerically solving for 

T. 

state solution very quickly (less than 1 second). 

Table [T] has an implication for modeling diffusion in porous media. At the pore scale the distances 
along pore throats ranges from microns to millimeters jlO| . This indicates that the water vapor along 
the pore can be viewed as in a quasi-steady state; that is, for any snapshot in time the relative humidity 
(similarly, concentration or chemical potential) is approximately given by the steady-state solution to (??). 
On the other hand, if we consider a porous medium where the gas phase is connected to atmospheric 
conditions, the diffusive length is increased by the tortuous nature of the connected gas-filled pores. Hence 
the transient equation may be used to model diffusion at this scale. The drawback (and the drawback 
to all pore-scale modeling) is that the tortuous nature of the diffusive path in the porous medium gives 
the mathematical model possibly intractable geometry. For this reason, upscaling techniques are needed. 
Through the upscaling, the diffusion coefficient becomes a function of the tortuosity, saturation, temperature, 
and possibly other variables. In this situation, one must not assume that the transient nature of the diffusion 
can be neglected. 



3.3 Sensitivity to Thermal Gradients 

To conclude the discussion on the diffusion dominated models we consider the sensitivity to thermal gradients. 
In 1957, Philip and deVries 8][9 18 suggested a possible model for enhanced vapor diffusion in porous media. 
This model states that Fickian diffusion is enhanced by pore-scale gradients in saturation and pore-scale 



gradients in temperature. The main argument against this model (eg 23 27]) is that pore-scale gradients are 
inherently immeasurable. Even so, we will show in this section that the relative humidity model is sensitive 
to thermal gradients. The introduction of a thermal gradient can drive diffusion faster than predicted by 
pure Fickian diffusion or can dampen diffusion so that the steady state solutions expected from Fickian 
diffusion do not hold. Throughout this subsection we assume that the gas is in thermal equilibrium with the 
liquid at the liquid-gas interface and that the temperature profiles are at steady state. 
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The saturated vapor pressure, Psat, is typically taken to only be a function of temperature. One empirical 
equation for the relationship between (absolute) temperature and saturated vapor pressure is the so called 
Antoine equation: 



exp (77.345 + 0.0057r + 7235/T) 

Psatyt) Jig 2 ■ 



(3.16) 



If we introduce a thermal gradient along the vertical axis of the tube, 



T(^,T) = (ri-To)C + To 



(3.17) 



then psat — Psat/iR^^T) is clearly a function of space and can no longer be considered constant. 

The diffusion coefficient, D, is also sensitive to temperature. According to Cussler jtj (equation 5.1-1 
page 104), the diffusion coefhcient is approximately proportional to T^^'^/p^: 



D cx 



7-3/2 



py 



(3.18) 



The proportionality constant is empirically fit and is a function of the types of molecules in the gas mixture. 
Using equations (??) and (??) the advection diffusion equation becomes 



dt 



V • ipsatiT)ip^n = V • [D{T)V {psatiT)ip)] 



(3.19) 



Let Dq be the diffusion coefficient at x = 0. Also let x = and t = [L'^ / Dq)t define the new dimensionless 
spatial and temporal groups respectively. With this, the dimensionless form of the advection diffusion 
equation is 



dT 



DjT) 

Dn 



{psat(T)ip) 



(3.20) 



As an experiment to determine the sensitivity of this equation to gradients in temperature, we induce 
a thermal gradient from the bottom to the top of the capillary tube (as in equation (??)) assuming a one 
dimensional geometry and Pe <^ 1. This closely mimics a pore-scale temperature gradient which deVries 18 
predicts as a mechanism for enhanced vapor diffusion in porous media. A baseline temperature of 300K is 
chosen with gradients of up to lOK. The solution to the PDE is found using MATLAB's pdepe routine. 

Figure |4] shows the relative L2 errors between the thermally forced diffusion and the linear diffusion model 
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as predicted with constant temperature 

Error [t; To, Ti) ^ — . 3.21) 

\mtr;To)\\2 

The boundary conditions were held fixed at 939 ~ 1 find ipi w 0.1 for this experiment as this represents 
an extreme case experimentally. For thermally forced diffusion we let Tq > Ti, and for thermally damped 
diffusion we let Tq < Ti. The names here are chosen to match the solutions seen in Figure [5) 



Relative Error: Thermally Damped vs. Linear (T^ = 300K) Relative Error: Ttiermally Forced vs. Linear (T^^ = 300K) 




10"^ 10"' 10° 10"^ 10"' 10° 

T = (D/L^)t 1 = (D/L^)t 



Figure 4: Relative errors for thermally forced diffusion (equation (??)) vs. linear diffusion (equation (??)) 
with Pe < 1. 

The sample solutions shown in Figure [5] illustrate the significant forcing (or damping) that a thermal 
gradient of lOK can cause. It is clear from Figure |4] that the thermally forced model predicts up to a 10% 
relative error against the linear diffusion model at steady state. This gives merit to the propositions by 
deVries in [Ts], but a pore-scale thermal gradient of lOK is unreasonable for most scales. A gradient of 
~ IK, on the other hand is possible in real media. For this reason we focus on the lower curves in Figure 
[4] In this case we see relative errors of approximately 1% for a one degree difference in temperatures. Table 
[2] shows the relative errors between the steady states as the gradient approaches zero. We see the expected 
convergence toward the linear diffusion model as Tq — ^ Ti or Ti Tq, but we also see a small enhancement 
even for very small gradient in temperature. 

The argument against deVries' model (pore-scale thermal gradients are inherently immeasurable) is still 
valid, but the present model also predicts thermally enhanced diffusion on a larger scale. If the porous 
medium has a thermal gradient across the entire diffusive (eg evaporative) length then the present model 
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Thermally Damped Diffusion: T =290K, T =300K 




Thermally Force Diffusion: Tg=300K, T^=290K 




Figure 5: Sample solutions of thermally damped (left) and forced (right) diffusion 





Ti = 299.0 


Ti = 299.5 


Ti = 299.9 


Ti = 299.99 


Ti = 299.999 


To ^ SOOK 


1.5% 


0.78% 


0.2% 


0.02% 


0.002% 




To = 299.0 


To = 299.5 


To = 299.9 


To = 299.99 


To = 299.999 


Ti = 300if 


1.5% 


0.77% 


0.2% 


0.02% 


0.002% 



Table 2: Relative L2 errors between thermally enhanced and linear diffusion steady states as VT — >■ 



will predict an enhancement. It should be noted that deVries also predicted enhanced vapor diffusion based 
on gradients in saturation. This problem is not handled in this work, but has been recently studied in 23125 



4 Advection Diffusion Models 



To conclude this work we now turn our attention to situations where Pe > 1. This physically corresponds to 
forced advection models. The prototypical example of a physical system governed by the advection diffusion 
equation is the diffusion of dye in a stream in the presence of a current. Relative to the present work this 
would be the diffusion of the water vapor through the capillary tube in the presence of wind (or some other 
forced convective current). This wind can have the effect of either forcing or damping the diffusion within 
the tube. 

To model this behavior we couple the advection diffusion equation, (??), with a Navier-Stokes type 
equation to govern the bulk gas-phase velocity. The focus of this work is to better understand the diffusion 
problems, so the equations and systems in this section will be stated without solution. 

To derive the proper form of the Navier-Stokes equation let us start with the momentum balance equation 
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for the gas 

ps^+p^vs- Vvs- V-f -p5g = 0. (4.1) 

Here is the Cauchy stress tensor for the gas phase. The Stokes assumption says that the stress tensor is 
the sum of the static pressure and viscous forces acting on the symmetric part of the gradient of velocity. 

r = -p^I + MVvf,„. (4.2) 

With Stokes' assumption we arrive at the traditional Navier-Stokes equation 

_ + p9^g . vv« + Vp5 - iiV ■ ( Vvf^„) - p»g = 0. (4.3) 

In the present situation there are two main cases to consider: (1) if Pe « 1 then we may assume that the 
viscous forces are negligible and the stress tensor can be approximated by the pressure and (2) if Pe » 1 
then forced convection is occurring and the viscous forces are likely non-negligible. Case (2) yields the full 
Navier-Stokes equations as seen in (??). In this work wc focus on natural convection so case (1) is the only 
relevant case. Consider the initial condition for to be zero throughout the domain and consider a no flux 
boundary condition along the sides of the capillary tube and along the liquid-gas interface. It remains to 
relate the pressure term to the relative humidity. 

In the absence of viscous forces, the Cauchy stress tensor is approximated by « — P^I- In this case the 
momentum balance equation, (??), becomes 

^ + p^vs • Vvs + VpS - pSg = 0. (4.4) 

We first consider a variable gas-phase pressure and then will simplify considering constant pressure. The 
assumption of constant pressure is the most restrictive, but on very small scales it may be reasonable. 

Under the assumption that the gas phase pressure is variable the ideal gas law can be used to rewrite the 

pressure in terms of the relative humidity. The pressure in the gas phase is the sum of the partial pressures 
of the water vapor and air (Dalton's Law), and the partial pressure of the water vapor can be rewritten in 
terms of the relative humidity. This is simply a statement of the ideal gas law: 

p9„ = ps^TpS" = RS^Tp,at{T)ip. (4.5) 
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Similarly, the partial pressure of the inert air can be written as 



pB'^=R3-Tff\ (4.6) 

but converting directly from air-species density to relative humidity is not straight forward. We choose to 
use the mass balance equation for the air-species to help close the system: 

_^+V.(p5«v«<')=0. (4.7) 

This introduces another closure issue since now an equation for the air-species velocity, v^", is needed. The 
equation balancing the upward motion of the water molecules (in an evaporative process) and the downward 
motion of the air molecules is given by equation (??). Multiplying by p^, using Pick's Law to replace p9vff9v,9 ^ 
and solving for pS'^v"" gives 

p9a^ga = DVp^^ + pS"v9. (4.8) 

Plugging (??) into (??) (and assuming that the temperature is constant) closes the system of equations for 
the case where is variable: 

^ + V • ((^vS) = DV -Vifi (4.9a) 

^ + V (p5»v«) = -Dps^tV ■ (4.9b) 

^ + pSv^ • Vv9 + VpS = p»g (4.9c) 

= p9v + = jiQv rpp^^^ (T)ip + Tps» (4.9d) 

= Psat{T)<fi + p'''. (4.9e) 

Equations (??) and (??) are mass balance equations for the relative humidity and the air-species mass 
density respectively. Equation (??) is a simplified form of the Navier-Stokes equation for the bulk gas 
velocity. Equation (??) is Dalton's Law law with the ideal gas law, and equation (??) is derived from the 
fact that the sum of the species densities is the bulk density. 

To write this system in dimensionless terms we must consider the particular system that we're trying 
to model. In typical analysis of the Navier-Stokes equations, an experimentalist is interested the Reynolds 
Number; the ratio of inertial forces to viscous forces. In the present case we have neglected viscous forces 
already and are more interested in how natural convection occurs as compared to diffusion. For this reason 
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we nondimensionalize the mass transport equations in the same manner as before: 



dr 



dp3 



(4.10a) 
(4.10b) 



Here we have defined p^" as the dimensionless air-species density via = psat{T)p^'' . Using the same 
dimensionless time and length as in equations (??) and (??), equations (??) - (??) become 











a7 + 





(4.11) 



In the present case we assume that Pe w 1 so from the definition of the Peclet mmiber (equation (??)), 
Vc ~ DjL. Therefore, 



((^ + 



1?) Ih 



L 



L 



V^^ - ( ^ ) V^pS" 



(4.12) 



Dividing by {W^T)/L (since we are most interested in the behavior of ip) we arrive at the dimensionless 
form of (??) - (??): 



L^Rs-T dr \L'^R9-T 



V9 . VfV^ - ( I 63 

« \Rg.T ' ^ 



R9- 
R9^ 



7fp9^. (4.13) 



For water vapor in a binary gas mixture, the diffusivity, D, is on the order of 10 ^ [m^/sec] [?]. Therefore 
the coefficient of the velocity terms is 



^L^R9-T 

and the coefficient of the gravitational term is 

f Lg 



0(10"^^), 



(4.14) 



\R9-T 



LO{10-^). 



(4.15) 



The ratio of the specific gas constants is 



R9- 



286.9 
461.5 



0.622. 



(4.16) 



Given the orders of magnitude of the dimensionless groups it is clear that the velocity terms are likely 
negligible for scales down to a hundredth of a micron and the gravity terms are likely always negligible. 
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Thus the momentum equation can be seen as an expression for the gradients of Lp and 



R9- 



0. 



(4.17) 



This is a statement about natural convection in a diffusive process: motion of the bulk gas phase is only due 
to gradients in density. Simply put, coupling a naturally convective system with the momentum equation 
yields a negligibly small bulk velocity. This implies that equation (??) is the only important equation and 
in the absence of forced convection we can safely solve the linear diffusion equation without regard to the 
momentum equation. 

In the case where is approximately constant throughout the domain, we can simplify and nondimcn- 



sionalize system (4.9) to 



d(p 



(Pe)V^ • {(pw^ 



v5 . Vfvf = 



63- 



(4.18a) 
(4.18b) 



Given that D ^ 0{1Q^^) we see that the graviatational term domanates the velocity terms for most length 
scales. This implies the same conclusion as in the case of variable pressure; without significant forcing, the 
bulk velocity term is likely negligible. Natural convection currents within the bulk gas are described through 
equation (??). 



5 Conclusions 

When considering diffusion in a simple capillary tube geometry there are three natural dependent variables in 
which to frame the problem: mass density, relative humidity, and chemical potential. The chemical potential 
has a natural description as a driving force for diffusion, but the transport equation only solves for it indirectly. 
A natural choice for modeling is to use the relative humidity as the dependent variable. In converting to 
relative humidity one must take care of the fact that the scaling factor, psat, is a function of temperature. 
Thermal gradients along the length of the capillary tube present an enhancement (or dampening) effect on 
the diffusion, and equation (??) can be seen as an alternative description of the enhancement factor predicted 
by deVries. Gradients in temperature on the order of IK over the length of the tube yield enhancements of 
roughly 1% L2 relative error in the transient behavior of the solutions as compared to the linear diffusion 
equation without thermal forcing. 

In this work, we found that a classical diffusion model by Bird, Lightfoot, and Stewart [4] may be mis- 
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interpreted so as to give enhanced vapor diffusion. This enhancement is a manifestation of neglecting the 
velocity of the air species within the binary mixture. Under the assumption that the temperature and gas 



phase pressure are approximately constant, we showed that the equations presented in 4|29 can be simplified 
back to the linear diffusion equation. 

When considering full advection diffusion equations one can still model in terms of the relative humidity 
(or equivalently, the mass conentration or chemical potential). Care must be taken in these models to pick 
a proper form of the Navier-Stokes equation (momentum balance) to model the system of interest. If the 
convective velocity is on the same order as the diffusive velocity (i.e. Pe ~ 1), then a simplified form of the 
Navier-Stokes equation can be used that neglects viscous forces. This coupling simply yields a description 
of natural convection current; hence indicating that the bulk velocity term is negligibly small in the absence 
of some significant forcing. 

A curious result of this work, as stated above, is that the chemical potential does not yield a more 
mathematically natural set of transport equations for diffusion dominated systems. This is contrary to the 
typically accepted thermodynamic definition of chemical potential as a descriptor for molecular diffusion. 
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A Nomenclature For Pore-Scale Advection Diffusion Models 

Superscripts, Subscripts, and Other Notations 





component of a— phase 




a— phase 




: difference of two quantities, = (•)" — (•)'' 


• a: 


(bold symbol) vector quantity 


• A: 


(double underline) second order tensor 


• (•)- 


a reference quantity or a quantity evaluated at a reference state 


• Oo: 


a quantity evaluated at the liquid-gas interface 


• Oi: 


a quantity evalued at the top of the capillary tube 


• (•): 


a dimensionless quantity 


Latin Symbols 


• A: 


dimensionless group A = (Ai?^" ) / i?^" [-] 




Mass fraction of j^^ compontnent C"-' = p"^ / p" [-] 


• D: 


diffusion coefficient [m^/s] 


• g: 


gravity [m/s^] 


• 33^: 


flux of species j [kg/ (m^-s)] 


• L: 


characteristic length (capillary tube length) [m] 


• Pe: 


Peclet number Pe = {Lvc)/D [-] 


• 


partial pressure of species j in phase a [N/m^] 


• p": 


pressure of a phase [N/m^] 


• Psat- 


saturated vapor pressure [N/m^] 


• R3i: 


specific gas constant for species j in the gas phase [J/(kg-K)] 


• T: 


absolute temperature [K] 
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• t: time [s] 

• t": Caucy stress tensor for the a phase [N/m^] 

• u: dimensionless chemical potential u = (/i^" — iJ,i")/{R^''T) [-] 

• v"^ : velocity of species j in phase a [m/s] 

• v": velocity of phase a [m/s] 

• v^: dimensionless gas phase velocity = Vc'v^ [-] 

• Vc- characteristic advective velocity [m/s] 

• x: geometric length [m] 
Greek Symbols 

• k: mean curvature of liquid-gas interface [1/m] 

• A: dimensionless group A = Psat/p» [-] 

• V|: gradient operator in terms of dimensionless length [-] 

• /u: dynamic viscosity [Pa-s] 

• fi^": chemical potential of water vapor in gas phase [J/kg] 

• /uf": chemical potential of water vapor at standard temperature and pressure [J/kg] 

• (fi: relative humidity (divided by 100) = p^'' /Psat [-] 

• p^i : mass density of species j in phase a [kg/ m^] 

• p": mass density of phase a [kg/m^] 

• psat' mass density of water vapor under saturated conditions [kg/m^] 

• a: surface tension [mN/m] 

• r: dimensionless time t = {Li^ /D)t [-] 

• 5: dimensionless length x = [-] 
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